Enrichment of the r-process Element Europium in the Galactic Halo 
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ABSTRACT 

We investigate the enrichment of europium, as a representative of r-process elements, 
in the Galactic halo. In present chemical evolution models, stars are assumed to be 
formed through shock processes by supernovae (SNe). The enrichment of the interstellar 
medium is calculated by a one- zone approach. The observed large dispersions in [Eu/Fe] 
for halo stars, converging with increasing metallicity, can be explained with our models. 
In addition, the mass range of SNe for the r-process site is constrained to be either stars 
of 8 - 10M Q or £ 3OM . 

Subject headings: Galaxy: evolution — Galaxy: halo — nuclear reactions, nucleosyn- 
thesis, abundances — stars: abundances — stars: Population II — supernovae: general 
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1. Introduction 

Recent abundance analysis of metal-poor halo stars 
reveals the presence of large dispersions in heavy ele- 
ments. This may be interpreted as a result of incom- 
plete mixing of the interst ellar medium (ISM ) at the 



beginning of the Galaxy (Gilroy et al. 1988). Each 



type of element shows a unique dispersion, which can- 
not be simply explained by spatial inhomogeneity of 
the ISM. The dispersions of neutron-capture elements 
like Sr, Ba, and Eu range ~ 300-fold, while those of a- 
elements and iron-peak eleme nts range typically <; 10 



fold (Mc William et al. 1995; Ryan, Norris, & Beer; 



1996). This fact implies mixing of ejecta from small 



numbers of SNe into the parent clouds (Audouze & 



Bilk 1995| ; Mc William 1997, 1998). 

In previous chemical evolution models, observed 
stellar compositions are taken to represent those of 
the ISM when the stars were formed. It may not be 
true, however, if star formations are mainly triggered 
by SNe. The composition of the formed star must be 
a mixture of the ISM and the individual SNe ejecta. 
Mathews, Bazan, & Cowan (1992) have examined the 
enrichment of neutron-capture elements using a chem- 
ical evolution model of the ISM. They concluded that 
the delayed increasing of [Eu/Fe] in terms of [Fe/H] 
favored origins of low mass SNe (~ 7 — 8M Q ). How- 
ever, this might change when including the large dis- 
persion in [Eu/Fe] for halo stars. 

The excellent agreement of the neutron-capture el- 
ements in the halo stars CS 22 892-052, HP 115444 
HD 122563, and HD 126238 flSncdcn et al. 1996; , 
1998) with the solar r-process abundance patt ern im- 
plies the pre sence of one robust r-process site flCowan 



]pt al. 1998 ). However, the origins of r-process el- 
ements are still unknown. Although the neutrino 
winds in SNe have been thought to be a promis- 
ing site, this scenario involves serious problems, e.g., 



in obtaining sufficient entropy ( 


Takahashi, Witti, & 


lanka 1994; 


Woosley et al. 1994; 


Qian & Woosley 



1996 ). Collapsing O- Ne-Mg cores are also thought to 
be an r-process site (Wheeler, Cowan, & Hillcbrandt 
1998). The O-Ne-Mg core, as the final stage of an 



8 — 10M Q star, may explode by a prompt shock (Hillc- 
^randt, Nomoto, fc Wolff 1984j ). The ejected shell 
contains a rather low electron fraction region due to 
electron captures, which may be a promising r-process 
site ( [Meyer fc Brown 1997| ). 

In this Letter, we discuss the enrichment of eu- 
ropium in halo stars formed through stimulated pro- 



cesses by SNe. Details of our models including other 
elements will be presented in a forthcoming paper. 

2. Methods 

The evolutions of the ISM in the Galactic halo are 
calculated by a one-zone halo model which loses gas 
through accretion onto the disk. The star formation 
and accretion rates are assumed to be proportional 
to the gas fraction of the halo. The star formations 
obey the Salpeter initial mass function in the mass 
range 0.05 — 60A/©. The coefficients for the accretion 
rate and the star formation rate are adjusted to fit to 
the observational data of [O/Fc] vers us [Fe/H] (e.g., 



Barbuy 198S ; Edvardsson et al. 1993) and the mctal- 
licity distribution of halo stars (Sandage & Foutg 



1987). Stellar lifetimes are adopted from Schaller et 
al. (1992). Type la SNe are supposed to occur simply 
with a lifetime of 2.5 Gyr for ~ 5% of the 3 - 8M Q 



1996) 



stars (ishimaru 1994; Yoshii, Tsujimoto, & Nomoto 



We assume that star formation is initiated by SNe. 
An SN remnant is supposed to expand spherically un- 
til reaching the merge radius with the ISM (typically 
- 100 pc; pioffi, McKee, fe Bertschinger 1988] ). The 
chemical composition of a formed star is assumed to 
be the mass average of the "snowplowed" ISM and 
the SN ejecta. The mass of the SN progenitor is cho- 
sen randomly but obeying the initial mass function. 
Further, we presume the following two cases in which 
the merge radii have a Gaussian distribution in the 
logarithmic scale within a factor of 1.5 (case 1) and 
no distribution (case 2). 

Stellar yields for Type II and Type la SNe are 
taken from Nomoto et al. (1997a) and Nomoto et 
al. (1997b), respectively. The 8 — 10M© stars are as- 
sumed to produce no iron. For simplicity, we do not 
take into account the metallicity effects as examined 
by Woosley and Weaver (1995). We do not include 
yields of the stars smaller than 8Mq either, which 
does not effect the results of the present study. The 
r-process elements are supposed to be produced only 
in Type II SNe. In this study, we examine the fol- 
lowing three models (hereafter models 1, 2, and 3) 
in which europium is produced from the stars: (1) 
8 - 10M© , (2) > 10M©, and (3) > 3OAf . For each 
model, the mass of produced europium is assumed to 
be constant over the stellar mass range, scaled to be 
[Eu/Fe] = 0.5 in the ISM at [Fe/H] = -2. 
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3. Enrichment of europium in the halo 

Figure 1 shows observational data of [Eu/Fe] versus 
[Fe/H] for halo and disk stars. Figures 2, 3, and 4 
illustrate the enrichment of europium in the halo by 
models 1,2, and 3, respectively. The evolution of the 
ISM is represented by a line, and the compositions of 
stars are plotted by open stars (case 1) and filled stars 
(case 2). As can be seen in the decrease of [Eu/Fe] for 
[Fe/H] J> — 1 in Figures 2-4, the iron production by 
Type la SNe has little contribution to the enrichment 
for metal-poor halo stars. 

We see from models 1 and 3 (Figs. 2 and 4) that 
a large dispersion in [Eu/Fe] appears for stars at 
[Fe/H] ~ — 3 : converging towards the value of the 
ISM with metallicity. These results are in excellent 
agreement with the observations (Fig. 1). The over- 
production of europium in stars compared to the value 
of the ISM is due to the star formations by SNe that 
add europium with little or no iron to the ISM. In 
contrast, the underproduction is due to SNe that add 
iron without europium. Thus, for model 1, a signifi- 
cant overabundance in [Eu/Fe] appears for the stars 
formed by 8-10M Q SNe. For model 3, stars > 30M© 
are supposed to produce europium with significantly 
high ratios relative to iron, since these stars account 
for only ~ 15% of all SNe in number. As a result, a 
large dispersion also appears. In contrast to models 1 
and 3, stars in model 2 (Fig. 3) show little dispersion 
in [Eu/Fe], since all SNe (> 1OM ) produce iron and 
europium with similar ratios. This result is rather 



[Fe/H] 

Fig. 2. — Enrichment of europium as a function of metallic- 
ity in the Galactic halo for model 1. The production site of 
europium is assumed to be 8 — IOMq stars. The evolution of 
the ISM is indicated by a line, and the compositions of stars 
are denoted by open stars (case 1) and filled stars (case 2). 
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Fig. 3. — Same as Fig. 2, but the site of europium is assumed 
to be > 1OM stars (model 2). 



close to the observational data of a-elements. 

The distribution of merge radii of SN remnants 
(case 1) can be another reason for the dispersion in 
[Eu/Fe]. For case 2, stars show a smaller scatter due 
to no distribution in the merge radii, and are thus 
clearly separated into overabundant and underabun- 
dant ones. 

As can be seen in Figures 2 and 4, the number of 
stars overabundant in europium is smaller than those 
that are underabundant. This is also in good agree- 
ment with the observational results (Fig. 1). The rea- 
son is simply that the number of SNe producing eu- 
ropium is smaller than of other SNe in these models. 



stars (Case 1 ) A 
stars (Case 2) * 
gas — 



[Fe/H] 

Fig. 4. — Same as Fig. 2, but the site of europium is assumed 
to be > 30Mq stars (model 3). 



4. Discussion and Conclusions 

In this study we find that a large dispersion in 
[Eu/Fe] for halo stars is mainly due to star formation 
by individual SNe that eject unique yields. There may 
be other reasons for the dispersion (e.g., a distribution 
of merge radii of SN remnants with the ISM [see §3] 
and a spatial inhomogeneity of the ISM), which are 
not dealt with in this study. We should emphasize, 
however, that the stochastic star formation process 
triggered by different masses of SNe is essential to re- 
produce the ~ 300-fold dispersions in neutron-capture 
elements. A dispersion caused by other factors may 
be ~ 10-fold at most, since it must appear similarly 
in other types (e.g., a- and iron-peak) of elements. 

The results of our study indicate that the produc- 
tion sites of europium, as a representative of r-process 
elements, must satisfy at least one of the following 
two conditions: (1) europium is produced with little 
iron or (2) the number of SNe producing europium ac- 
counts for only a small fraction of all SNe. We suggest 
here a couple of possible sites for these conditions. 
The first is the explosion of 8 — 1OM stars. These 
stars are expected to produce little iron (Hillebrandt 
et al. 1984). In fact, the progenitor of the Crab Neb- 
ula, which shows no significant enrichment in metal, 
ha s been su spected to be an 8 — 1OM star ( Nomoto 
et i J. 1982| ) . The second site is the explosion of stars 
;> 30-Mq. The subsequent neutrino winds may ob- 
tain substantially high entropy owing to the massive 
neutron stars ( Qian fc Woosley 1996 ). These stars 
account for only ~ 15% of all SNe. Furthermore, the 
explosions of these stars may eject significantly less 




iron than ~ O.IMq adopted in this study (Woosley & 
Weaver 1995). In fact, very low masses of 56 Ni in the 
ejecta of SN 1994W (~ 0.0026Af M ; |Sollerman, Cum- 



ming, fc Lundqvist 1998| ) and SN 1997D (~ 0.002M Q ; 

1998) have been expected to be due to 



25 — 40Mq progenitors. It should be noted that coa- 
lescing neutron stars may also produce r-process ele- 
ments without iron ( Rosswog et al. 1998| ). However, 
they may not be major sources for the large enhance- 
ment in [Eu/Fe], since their lower kinetic energies are 
not enough to trigger star formations. 

In the light of the observational results (Fig. 1), 
model 1 seems more likely for the following two rea- 
sons. First, in model 1 some stars show as high an 
[Eu/Fe] as ~ 1.7 at [Fc/H] ~ -3. These stars can 
be a possible explanation for CS 22892-052. How- 
ever, model 3 would produce stars with a significantly 
larger dispersion in [Eu/Fe] if smaller iron yields were 
adopted for stars > 30M©. Second, in model 1 
there is no sharp peak in the distribution of stars 
for [Fe/H] < — 2 (case 1) that is in good agreement 
with the observations. In model 3 the concentration 
of stars on the line of the ISM for [Fe/H] < -2 is 
due to the star formations by 8 — 10M Q stars which 
add neither iron nor europium. It should be noted, 
however, that the question whether 8 — 10M Q stars 
explode or not is still open. Another outstanding dif- 
ference between these models is that model 1 has stars 
in [Eu/Fe] down to ~ — 2 that have not detected by 
observations yet. This is a consequence of the delayed 
increasing of europium in the ISM due to lower mass 
progenitors. Hence, future surveys of metal-poor halo 
stars for low [Eu/Fe] will be important to distinguish 
the above two sites for the r-process. An alternative 
way to distinguish these sites may be to apply our 
models to other neutron-capture elements like Sr and 
Ba which have much more data than Eu. In addition, 
the decreasing trend of [Ba/Fe] towards lower metal- 
licity may support 8 — 10M Q progenitors. The future 
direction of this study will be one that investigates 
enrichment of a number of neutron-capture elements 
consistently. 

We would like to acknowledge useful discussions 
with T. Young, S. Ryan, and T. Kajino. This work 
was supported in part by Japan Society for Promotion 
of Science. 
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